Charged meson rapidity distributions in central Au+Au collisions at y / J/\F/v = 200 GeV 
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We have measured rapidity densities dN/dy of tt^ and K ± over a broad rapidity range (—0.1 < 
y < 3.5) for central Au+Au collisions at ^/snn = 200 GeV. These data have significant implications 
for the chemistry and dynamics of the dense system that is initially created in the collisions. The full 
phase-space yields are 1660±15±133 (tt+), 1683±16±135 (tt _ ), 286±5±23 (K+) and 242±4±19 
{K~). The systematics of the strange to non-strange meson ratios are found to track the variation of 
the baryo-chemical potential with rapidity and energy. Landau-Carruthers hydrodynamic is found 
to describe the bulk transport of the pions in the longitudinal direction. 



In ultra-relativistic heavy ion collisions at RHIC en- 
ergies, charged pions and kaons are produced copiously. 
The yields of these light mesons are indicators of the en- 
tropy and strangeness created in the reactions, sensitive 
observables to the possible existence of an early color 
deconfined phase, the so-called quark gluon plasma. In 
such collisions, the large number of produced particles 
and their subsequent reinteractions, either at the par- 
tonic or hadronic level, motivate the application of con- 
cepts of gas or fluid dynamics in their interpretation. Hy- 
dro dynamical properties of the expanding matter created 
in heavy ion reactions have been discussed by Landau Q 
(full stopping) and Bjorken |2j (transparency), in theo- 
retical pictures using different initial conditions. In both 
scenarios, thermal equilibrium is quickly achieved and 
the subsequent isentropic expansion is governed by hy- 
drodynamics. The relative abundances and kinematic 
properties of particles provide an important tool for test- 
ing whether equilibrium occurs in the course of the colli- 
sion. In discussing the source characteristics, it is impor- 
tant to measure most of the produced particles in order 
not to violate conservation laws (e.g. strangeness and 
charge conservation). 

In this letter, we report on the first measure- 
ments at RHIC energies of transverse momentum 



(pt) spectra of n ± and if ± over the rapidity range 
—0.1 < y < 3.5 for the 5% most central Au+Au collisions 
at y/sjvN = 200 GeV. The spectra are integrated to 
obtain yields as a function of rapidity (dN/dy), giving 
full phase-space (47r) yields. At RHIC energies, a low 
net-baryon density is observed at mid-rapidity 
so mesons may be predominantly produced from the 
decay of the strong color field created initially. At 
forward rapidities, where primordial baryons are more 
abundant other production mechanisms, for example 
associated strangeness production, play a larger role. 
Therefore, the observed rapidity distributions provide 
a sensitive test of models describing the space-time 
evolution of the reaction, such as Landau and Bjorken 
models 0, 0. In addition, integrated yields are a key 
input to statistical models of particle production [5J, |6j . 



BRAHMS consists of two hadron spectrometers, a 
mid-rapidity arm (MRS) and a forward rapidity arm 
(FS), as well as a set of detectors for global event char- 
acterization J7J. Collision centrality is determined from 
charged particle multiplicities, measured by scintillator 
tile and silicon multiplicity arrays located around the 
nominal interaction point. The interaction vertex is mea- 
sured with a resolution of 0.6 cm by arrays of Cerenkov 
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counters positioned on either side of the nominal vertex. 
Particle identification (PID) for momenta below 2 GeV/c 
is performed via time-of-flight (TOF) in the MRS. In 
the FS, TOF capabilities allow ir-K separation up to 
p = 4.5GeV/c, and is further extended up to 20GeV/c 
using a ring imaging Cerenkov detector. Further details 
can be found in 0, • 

Figure ^ shows transverse mass mx — mo spectra 
(rriT — \J V \ + m o) f° r 7r_ an d K~ . Particle spectra 
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FIG. 1: Invariant transverse mass mr — mo spectra of ir~ (a) 
and K~ (b) from y ~ (top) to y ~ 3.5 (bottom). Dashed 
lines are fits to the data, namely a power law in px for pi- 
ons and an exponential in mr — mo for kaons. Errors are 
statistical. Spectra have been rescaled by powers of 10 for 
clarity. 

were obtained by combining data from several spectrom- 
eter settings (magnetic field and angle), each of which 
covers a small region of the phase-space (y,px)- The 
data have been corrected for the limited acceptance of 
the spectrometers using a Monte-Carlo calculation simu- 
lating the geometry and tracking of the BRAHMS detec- 
tor system. Detector efficiency, multiple scattering and 
in-flight decay corrections have been estimated using the 
same technique. Hyperon (A) and neutral kaon Kq s de- 
cays may have contaminated the pion sample. For Kq s , 
it is assumed that its yields amount to the average be- 
tween K + and K~ at each rapidity interval. For A yields, 
since only mid-rapidity data are available, we used the 
same assumptions as in 0, namely A/p — A/p — 0.9 in 
the phase-space covered in this analysis. The fraction of 
pions originating from A and Kq s decays was estimated 
with a GEANT simulation where realistic particle distri- 
butions (following an exponential in mx) were generated 
for several spectrometer settings. Particles were tracked 
through the spectrometers and produced pions were ac- 
cepted according to the same data cuts applied to the 
experimental data. It has been found a Kq s (A) contam- 
ination of 4% (< 1%) in the MRS and 6% (< 1%) in the 



forward spectrometer. In the following, the pion yields 
are corrected unless stated otherwise. 

The pion spectra are well described at all rapidities by 
a power law in px, A(l + px / Po) n ■ For kaons, an ex- 
ponential in mr — mo, Aexp ( mT ~ m ° ), nas been used. 
The invariant yields dN/dy were calculated by integrat- 
ing the fit functions over the full pt or mj- range. The 
two main sources of systematic error on dN/dy and (px) 
are the extrapolation in the low px range outside the ac- 
ceptance, and the normalization of the spectra. Other 
fit functions were used in order to estimate the error on 
the extrapolation. In the FS, due to a smaller acceptance 
coverage at low px , the error is systematically larger than 
in the MRS. In total, the systematic error amounts to ~ 
10% in the range -0.1 < y < 1.4 (MRS) and - 15% 
for y > 2 (FS). Mid-rapidity yields recently reported by 
the STAR and PHENIX experiments [ijj are within 
1 Psyst of these results. 

Rapidity densities and mean transverse momenta (px) 
extracted from the fits are shown in Fig. |2 Panel (a) 
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FIG. 2: Pion and kaon rapidity densities (a) and their mean 
transverse momentum {px) (b) as a function of rapidity. Er- 
rors are statistical. The kaon yields were multiplied by 4 for 
clarity. The dashed lines in (a) are Gaussian fits to the dN/ dy 
distributions (see text). 

shows the pion and kaon yields. 7r + and ir~ are found in 
nearly equal amounts within the rapidity range covered, 
while an excess of K + over K~ is observed to increase 
with rapidity 0. FigureJ^b) shows the rapidity depen- 
dence of (px) ■ There is no significant difference between 
positive and negative particles of a given mass. For pi- 
ons, (p T ) = 0.45 ± 0.05 GeV/c (stat + syst) at y = and 
decreases little to 0.40 ± 0.06 GeV/c at y = 3.5, while for 
kaons, (px) drops from 0.71 ± 0.07GeV/c at y = to 
0.59 ± 0.09GeV/c at y = 3.3 (see [n|). 

In order to extract full phase space densities for ^ 
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and K we have investigated several fit functions: a sin- 
gle Gaussian centered at y = (Gf), a sum of two Gaus- 
sians (G2) or Woods-Saxon (WS) distributions placed 
symmetrically around y = 0. Our data do not distin- 
guish among these functions, all give a x 2 per d.o.f. of 
~ 1, and have the same total integral to within 2% (cf. 
Tab.0. 
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TABLE I: Full phase-space yields of and K extracted 
from fits to the dN/dy distributions (see text). Errors are 
statistical, systematic errors are of the order of 8%. 

Figure |3f a) shows ratios of the strange to non-strange 
full phase-space meson yields (K/tt) as a function of 
^snn- The ratio K + /ir + shows a fast increase from low 
AGS to low SPS energies (^/snn = 8 GeV), followed by 
a decrease with increasing energy. At ^/snn — 200 GeV, 
we find a value of 0.173 ± 0.003 ± 0.014, consistent 
with the ratio at the highest SPS energy. In contrast, 
K~/t:~ increases monotonically but remains below 
K + /ir + . At y?SNN — 200 GeV, it reaches a value of 
0.143 ± 0.003 ± 0.011, a value close to the positive 
ratio. It is interesting to compare the energy systematic 
with the rapidity dependence at ^/snn = 200 GeV 
(Fig. EJb)). A fit to a straight line in the range 
< y < 1.3 gives K+ /tt+ = 0.162 ± 0.001 (stat) 
and K-j-K- = 0.154 ± 0.001. The dashed and dot- 
ted lines are predictions of the hadron gas statistical 
model 0, where a chemical freeze-out temperature T 
of 177 MeV and baryo-chemical potential /.ib of 29MeV 
were used (the authors used a fit to particle ratios 
at y/SNN — 130 GeV in the rapidity range \y\ < 0.5 
for extrapolating mid-rapidity ratios at the top RHIC 
energy). The agreement with the data is consistent 
within the systematic errors j2^. However, the data 
deviate from the model prediction at y > 2, where an 
increasing excess of K + over K~ is observed. This may 
be due to an increase of the net-baryon densities 0, 13 • 
A baryon rich environment is favorable for associated 
strangeness production, e.g. p + p — > p + K + + A, a 
production channel forbidden to K~ . In the context of 
the statistical model, this translates into an increase of 
the baryo-chemical potential /is, as already reported 
in Q, where a calculation by Becattini et al. [j| of 
K + 1 K~ vs p/p at constant T (and varying /i^) agrees 
well with the rapidity dependence of the experimental 
ratios. 

The observed rapidity distributions of pions exhibit a 
nearly Gaussian shape. Widths are found to be av+ — 
2.25 ± 0.02 (stat) and er^- = 2.29 ± 0.02 rapidity units. 
Similar overall features have already been observed in 
central Au+Au collisions at AGS fijl and Pb+Pb reac- 



tions at SPS [T^. This is reminiscent of the hydrody- 
namical expansion model proposed by Landau In 
the initial state, colliding nuclei are highly Lorentz con- 
tracted along the beam direction. Under assumptions 
of full stopping and isentropic expansion after the ini- 
tial compression phase (where thermal equilibrium is 
reached), the hydrodynamical equations (using the equa- 
tion of state of a relativistic gas of massless particles) 
lead to dN/dy distributions of Gaussian shape at freeze 
out 0|- The original Landau model was simplified and 
used by Carruthers et al. for the description of particle 
production in p+p collisions [T^. I20I l2l| . In their model, 
the width a of the pion distribution is simply formulated: 

a 1 = ln7 fceam = In (y r s/2m p ) (1) 

where m p is the proton mass. To extend this formula 
to heavy ion collisions, one might naively replace y/s by 
^/sWn- Figure a) shows the experimental dN/dy(ir) 
and Carruthers Gaussians at ^/snn = 200 GeV (us- 
ing Eq. ^ w hh the condition that the integrals of these 
Gaussians must be equal to the full-space yields ex- 
tracted from the data). Surprisingly, a discrepancy of 
only ~ 5% between the model and the data is observed 
(ccarrut. = 2.16). The inset in Fig. 21 shows the ratio 
Cdata/ccarrut. as a function of ^s NN . At all energies, 
the widths (Td&t& result from Gaussian fits to the pion 
distributions. While the difference between theory and 
measurements is of the order of 10% at most from AGS 
to RHIC energies, it is worth noting that the overall sys- 
tematic may not be trivial. The ratio at RHIC is ~ 15% 
higher than at SPS. It is furthermore interesting to note 
that in central Au+Au collisions at ^snn = 200 GeV, 
the original baryons lose 2 rapidity units on average from 
the initial value yb = 5.36 Not only is the degree 
of transparency significantly different between AGS and 
RHIC, but the relative rapidity loss (Sy) /yb is about half 
lower Q . 

On the basis of the Landau's original hydrodynamic, 
Bjorken 0] proposed a scenario in which yields of pro- 
duced particles would be boost-invariant within a region 
around mid-rapidity. In that approach, reactions are 
described as highly transparent leading to a vanishing 
net-baryon density around mid-rapidity and particle 
production from pair creation from the color field in the 
central zone. This would result in a flat distribution of 
particle yields around y — 0. As mentioned, collisions 
at RHIC are neither fully stopped nor fully transparent, 
although a significant degree of transparency is observed. 
Consequently the overall dN/dy distribution of pions is 
expected to consist of the sum of the particles produced 
in the boost invariant central zone and the particles 
produced by the excited fragments. The fact that the 
observed distributions are flatter at mid-rapidity and 
wider than those predicted by the Landau-Carruthcrs 
model might point in this direction. 

In summary, we have measured transverse momentum 
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FIG. 3: Full phase-space lf/7r ratios as a function of ^/snn (a) and rapidity systematics at ^/sjvjv = 200 GeV (b). The dashed 
and dotted lines in (b) are predictions of the statistical mode l . Errors are statistical and systematics in (a), onl y st atistical 
in (b). AGS data are from [UEU, SPS data from [3 [H Gj|. Data points at ^1^7 = 6.3 GeV and 7.6 GeV jH are 
preliminary. 
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FIG. 4: Comparison dN/dy(ji) and Landau's prediction at 
y/SNN = 200 GeV (a) and ratio <Tjv(-7r)/crcarrut. as a function 
of y/SNN (b). Errors are statistical. 



spectra and inclusive invariant yields of charged meson 
tt^~ and K^. The ratios of strange to non-strange 
mesons K/ir are well reproduced by the hadron gas 
statistical model that assumes strangeness equilibra- 
tion at mid-rapidity. The excess of K + over K~ yields 
at higher rapidities can be explained by the increasing 
baryo-chemical potential /is with rapidity. The widths 
of the pion rapidity distributions are in surprisingly 
good agreement with a hydrodynamic model based on 
the Landau expansion picture. 
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